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The phenomenon of "overlapping resonances interference-induced transparency" (ORIT) is introduced and 
studied in detail for the Sq — > S2/S1 photoexcitation of cold pyrazine (C4H4N2). In ORIT a molecule be- 
comes transparent at specific wavelengths due to interferences between envelopes of spectral lines displaying 
overlapping resonances. An example is the S2 O Si internal conversion in pyrazine where destructive in- 
terference between overlapping resonances causes the So — > S2/S1 light absorption to disappear at certain 
wavelengths. ORIT may be of practical importance in multi-component mixtures where it would allow 
for the selective excitation of some molecules in preference to others. Interference induced cross section 
enhancement is also shown. 

I. INTRODUCTION 

In "Electromagnetically Induced Transparency" (EIT) l|-|8(, one creates a photoabsorption transparency win- 
dow at certain frequencies by applying a ("coupling") laser operating at another set of frequencies. In this paper 
we investigate a related, though distinct, physical phenomenon, called "overlapping resonances interference-induced 
transparency" (ORIT), where the transparency occurs due to interference between material waves within a molecule. 
Though ORIT is known for small systems P, |l^| , it has not been investigated for polyatomic molecules where over- 
lapping resonances are far more ubiquitous. In the present paper we introduce the ORIT phenomenon for polyatomic 
molecules, and examine it in detail in the Sq — > S2/S1 photoabsorption of pyrazine (C4H4N2). 
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II. THE PYRAZINE So -»■ S 2 /Si PHOTOEXCITATION 



The ultrafast dynamics of pyrazine is of longstanding interest. For example, in recent work where the compu- 
tational tools used in this paper were develo ped, we studied pyrazine S2 <-> Si internal conversion, using a pre-excited 
superposition state in S2 as a starting point 



laser excitation 



15 
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14[ , or creating such superposition state in the So — » S 2 / Si ultrafast 
16 1 . In the course of these prior studies [12] a detailed understanding of the vibronic structure 



and intramolecular dynamics of pyrazine was obtained. 

We denote as intramolecular \n) "resonances" the vibrational states of the noninteracting S2 electronic state, 
and the projector onto this resonance manifold as Q = ^2 K \k)(k\. These resonances are coupled to the vibrational 



states of the Si electronic state, denoted as \f3), with P = J2r be ing 



the projector onto the manifold of these 



17 ], the eigenstates of the excited Si + S2 



states. When we ignore the coupling to the pyrazine triplet T states 
Hamiltonian, denoted as (7), take into account the Si O S2 couplings. We thus have that P + Q = I = ^2 ItKtI 
where / is the identity operator on the excited states manifold. Even in the neglect of rotations, the So — > S2/S1 



photoabsorption process becomes a formidable 24 modes vibrational problem. 

We now examine conditions under which the (weak held) photoabsorption cross section a(E) given as 



a{Ky) 



i<7i/%>r, 



(i) 



vanishes, where \g) is the ground vibrational state of So, /x is the transition-dipole operator, and w 7i9 = (E 7 — E g )/h 
is the excitation frequency. For E 7 above the 5*2 potential minimum, the full excited vibronic state (7) is composed 
of both the S2 and the Si vibrational states (remembering that the Si potential minimum lies below the S2 potential 



minimum). It was previously established 
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13. 



191 ] that (k|/i|<7), the dipole matrix elements for the pure So 
transitions, are an order of magnitude larger that ((3\/i\g), the analogous matrix elements for the pure So 
transitions. Hence, one can introduce the so-called "doorway states" approximation, according to which, 

(tMs) = h\(P + QMg) =$>l/3)0%ls> +£<7|k><«Hs> « $>|k)<kHs>- 

/3 K K 

This yields, for cr(iJ 7 ): 

2 

^_2, . 

^(7|k)(k|/i| 5 ) . 



S 2 
Si 

(2) 



(3) 



The cross section in Eq. ([3]) can be expanded as 



a(E y ) 



The first (positive) part of Eq. ((3]), 



a di ^(i? 7 )=^ 



|( 7 | K )( K Hg)| 2 ^E^ ias (^)' 



(4) 



(5) 



is the diagonal contribution to cr(£ , 7 ), composed of individual contributions <7j? lag (£ 7 ) from the \k) resonances, while 
a"^(E y ) = £ ^ M <fflMl«)<«Vl5>(«l7)(7l« / ) = E <«' rf ( £; 7) = E 2 °%?( E J (6) 

is the (positive or negative) interference contribution to a{E 1 ) : composed of pairwise interference terms cr™^, {E 1 ) , 
k ^ k' (crj t n ^°/ rf (-Ely) in Eq. (|6]) are assumed real). Each pairwise interference term <rj c n ^ c / rf (-E* 7 ) is proportional to the 
product (re|7)(7|re'), which is non-zero only if both and (7|k') are non-zero. When this happens, the \k) and \k') 

resonances are said to overlap. The overlap is due to the common contribution of the full state I7) (Ref. 
In other words, for the pairwise interference contribution <7™^°, rf (-E) between resonances \n) and to be non-zero at 
the particular energy E = E~ n \k) and \k') must overlap at this energy through one or more common states I7). 



The cr(Ej) cross section is closely related to the S2 excited state population Ps 2 (t) [16[, generated by CW pulse 
e(lj) « e a 5{ijJ — w 7)S ) that excites only one state I7). To obtain Ps 2 (t), consider the excited wave packet |^l/ a (t)) on Si 
and S2, produced by the one-photon transition from So in the framework of first-order perturbation theory. For the 
CW case \^ a (t)} takes the form: 



|*a(t)> = ^expHf? 7 t/R)|7)<7| 



£<«i 



(7) 



Taking into account Eq. ([2]), we see that the last sum of Eq. is fi\g) = \^ e ) - an excited wave packet on S^- Thus, 

|*«(t)> = (ie a /h)e^p{-iE 7 t/n)\ 7 )(^\g) = (ze Q /^exp(-zi? 7 t/ft)| 7 )(7|* e ), (8) 

a result proportional to the action of the single-state propagator exp(— iE^t/h)]^)^] on |\I/ e ). Noting that the popu- 
lation in the S2 is given as Ps 2 (t) = a (t)\Q\^ a (t)) , we can use Eq. ((8]) to obtain that 



Ps 2 = J|(* e |7)(7lQl7)(7l* e ) = §l(7l^ e )| 2 (7lQl7) = ^Kl^lf 



£k«|7>i j 



(9) 



That is, the CW-generated S2 population is time-independent, but |7)-dependent. 

The cr(E 7 ) cross section of Eq. ([1]) can also be written as the expectation value of the single-state projector 
I7X7I in the fi\g) = |* e ) state, 

c c c 

Comparing Eqs. ([9]) and (fTUf gives the following relationship between the CW-generated population at E 1 and the 
cross section a{E 1 ): 
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III. COMPUTATIONAL RESULTS FOR THE PYRAZINE S -> S2/S1 PHOTOEXCITATION CROSS 

SECTION 



A. Coarse-Graining of the Pyrazine Vibronic Structure 



Pyrazine has 24 vibrational degrees of freedom and approximately 10 10 I7) states 



12, 



20j at -2 eV above the 



5„ - 5, vertical e„e r6 v difference. For the cot— s to be _ffv feasible, we replace the exact states with 
a set of approximate "coarse-grained" states [12j . We thus divide the energy axis into 2000 bins I a of size A Q , center 
energy E a and density of states p a . The projector onto a coarse-grained state \a) is defined as 

l"X«l = (l/(p Q A Q )) E l7><7l) hcncc v / PaA Q |a)(a| v / PaA Q = 51 It) <tI ■ 

Thus, the coarse-grained state \a) effectively replaces all the I7) states in bin I a . Numerically, the weighted states 



I a) = \/ p a A a \a) and their overlaps with resonances \k) are available through our iterative solution method for pyrazine 
(described in detail in Ref. |12j). 



The above coarse-grained description of the 



S2 <-> Si internal conversion dynamics 



12 
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pyrazine vibronic structure is fully adequate for vibronic femtosecond 



16j . It involves a manifold of 76775 \a) coarse-grained vibronic states, 
spanning the energy range of 2 eV above the So — > Si vertical electronic transition. The density of |a) states is 
nonuniform with average vibronic energy separation of 2.6xl0~ 5 eV = 0.21 cm -1 , which happens to coincide with 
some rotational energy spacings. Thus, if one considers vibronic Sq — > S2 excitation by a CW laser, which is narrow 
enough to resolve particular \a) states separately, one also needs to account for the rotational transitions and consider 
the full ro-vibronic spectrum. However, as discussed in the Appendix, if one assumes that pyrazine is initially cold, 
then consideration of vibronic levels only is sufficient. 

To compute the cross section numerically, we thus replace (7|k) by (a\K) and E 7 by E a , throughout. When 
this is done, the cross section for the state \a) is given [using Eq. (JTJ|] as 

a(E a ) = \(a\n\g)f = ^^ {g \p\a)(a\^\ g ) = = 1a!^| (s |^)| 2; (i 2 ) 

c c c c 

where uj a , g = (E a — E g )/h. In terms of I7) states, this gives 

a(E a ) = ^^(* e |7)(7l* e > * E ^^(* e |7)(7l* e ) = E ^ 

7E-fa 76-Tq 7G-f Q 

Hence, <j(E a ) is approximately the cumulative sum of all the individual cross sections for all the I7) in this bin I a . 
This is a reasonable approximation as long as the I a bin size is small. 

As in Eq. ([3]), the coarse-grained cross section of Eq. (fT2")) can be decomposed into diagonal and interference 
contributions, 

o d ^{E a ) = E \{a\n){n\p\g)\ 2 E^^(£ a ), (14) 



5 

and 

= J2 <"\E a ) = J2 2 <^ f (E a ), (15) 

respectively (cJ. n ^ e / rf (-E a ) in Eq. (IT5l) are real). 

This result is analogous to that obtained earlier. That is, each pairwise interference term a 1 ^, (E a ) in the 
coarse-grained case is proportional to the product (n\a)(a\K'), which is non-zero only if both (k|«) and (cx\k') are 
non-zero, i.e., the resonances \k) and \k') overlap by means of the common coarse-grained state \a) jl^. Ilq]. Only 
under such circumstances cr™^°, rf (E a ) is non-zero. 



B. Computational Results for the So — > S2/S1 Photoexcitation Cross Section 



In this study, transition dipolc matrix elements (K\fj,\g) are approximated by the corresponding Franck-Condon 
factors [12I ]. One hundred and seventy six "brightest" |k) resonances, having the largest Franck-Condon factors, were 
taken into account. These bright \k) are numbered hereafter in order of ascending zero-order energy E®. 

The interference term in Eq. (|15l) can be positive or negative at any given E a , resulting in a constructive or 
destructive interference contribution to the photoabsorption cross section. The latter case is of particular interest, 
because it can lead to interference-induced transparency in the absorption profile. As an example, we consider the 
neighborhood of the particular resonance \k = 94). The diagonal term af &s {E a ) for this resonance is shown in Fig. 
[TJ together with the position of its zero-order energy, E® =94 5.1342 eV. The energy of the resonance maximum, 
E™J$ 4 « 5.1350 eV, differs from £° =94 by a value A K=94 « 8.0xl(T 4 eV (6.45 cm" 1 ), a resonance shift due to S1S2 
vibronic nonadiabatic coupling. 

For the particular energy E™g 4 , the contribution cr^g 4 (E™j*:g 4 ) from the particular resonance \n = 94) domi- 
nates the full diagonal cross section er dia s (E™J% 4 ) , with a^ 4 {E^ 4 ) / a dias {E^ 4 ) ~ 0.80 (see Fig. [2]). One expects 
the total cross section a(E^g 4 ) to be of the same order of magnitude, but in this particular case the total cross section 
a(E™*g 4 ) is actually much smaller than both er dias (E™*g 4 ) and cr^f 4 (I?£!f 4 ), leading to the local transparency. This 
transparency is due to the large negative (destructive) interference contribution cr mterf (£™f* 4 ). The result is shown in 
Fig. O Namely, the total cross section, composed of contributions from 176 bright \n) resonances, is 23.4 times smaller 
at E = E^g 4 than the contribution from the one resonance \k = 94), an effect due to overlapping \k) resonances. 

It is instructive to analyze the structure of <J dme (E™g 4 ) and cr mtorf (E™g 4 ). The former is composed of 176 
terms crf ae {E^ 4 ), shown in Fig. [3] One can see that, out of overall 176, besides n = 94 itself, k = 40, 70, 84, 98, 
102 and 114 contribute significantly to cr diag (£^4) . This is a generic feature of pyrazine. That is, for any bright \k) 
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FIG. 1: Cross section contribution cr^| 4 (S), provided by the resonance \k = 94). Position of E®= 94 is shown as a vertical 
dashed line and an arrow on the Energy axis. 
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FIG. 2: Cross section contribution <r^ las (i?) of resonance \n = 94), together with full diagonal cross section contribution 
cr dias (E), full interference cross section contribution cr mtcrf (_E), and full cross section cr(E). Position of E™=g 4 is shown as a 
vertical dash-dotted line and an arrow on the Energy axis. 



7 



only several contributions are found to be important in cr dlag (-E™ ax ). 
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FIG. 3: Individual cross section contributions af^iE) to cr diag (S) at E = E™J% 4 , k = 1-176. 

For k = 1-176, there are 176(176-l)/2 = 15400 different pairwisc interference contributions cr" 1 ^ e , rf (E) in 
<7 lnterf These interference contributions are presented in Fig. @] for E = E™^g 4 , using a one-dimensional cu- 
mulative index k + k' (k < k') for simplicity. The most prominent contributions to the interference part of the cross 
section are seen to arise from pairs of resonances with largest (by absolute value) contributions to the diagonal part 
(Fig. [3j. Namely, the most significant pairwise interference contributions, seen in Fig. |4l are the ones with k + k! 
= 134 {i.e. 40 + 94), 164 (70 + 94), 178 (84 + 94), 192 (94 + 98), 196 (94 + 102) and 208 (94 + 114). Thus, the 
overlap of K — 94 resonance with these several resonances is crucial for the destructive interference causing ORIT at 

77 pmax 

E - E K=g4 . 

In addition to ORIT, the interference term in Eq. (1151) can be positive, thus leading to constructive interference, 
amplifying the total cross section <r(E). One characteristic example of such interference is shown in Fig. [S]for the 
coarse-grained vibronic state \a = 66984), having E a — 5.704590176 eV, for which total cross section cr(E a ) is 10.7 
times larger than the diagonal cross section term a dlas (E a ); in other words, the interference term <r lntcrf (E a ) here is 
9.7 times larger than cr dia s(.E Q ). While in case of ORIT \a intcri (E a )\ < a di ^{E a ) due to non-negativity of total cross 
section a(E a ) = cr dlag (E a ) + cr lntcrf (E a ) , there is no such a limitation on the absolute value of constructive interference 
term with respect to diagonal term, as Fig. [5] shows. 
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FIG. 4: Pairwise interference cross section contributions 2a^, Ti (E) to a intcli (E) at E = £™=f 4 , k, «' = 1-176 (k < k'). 
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FIG. 5: Full diagonal cross section contribution cr dlag (i5), full interference cross section contribution cr mtcrf (£) and full cross 
section cr(E). Energy of the coarse-grained state \a = 66984), -E a =66984 = 5.704590 1 76 eV, is shown as a vertical dash-dotted 
line and an arrow on the Energy axis. 
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IV. So ->• Sa/Si PHOTOEXCITATION CROSS SECTION ORIT AND CW-GENERATED S 2 

POPULATION 



The example of ORIT, presented and discussed above, is provided in terms of microscopic cr{E) cross section. 
Physically, these effects should be observed in terms of full vibronic states I7). However, the S2 O Si coupling, which 
is responsible for the very dense and rich a(E) structure, makes the experimental photoabsorption spectrum (obtained 
with finite spectral resolution and in the presence of environment) diffuse and continuous. This is the case for all 
pyrazine So — > S2 photoabsorption spectra, presented in literature 21[, so that these available spectra do not allow 



comparison with our ORIT computation. So, instead, in this section we computationally confirm the relevance of 
ORIT in the cross section to the amount of CW-generated S2 population, while the comparison with the experimental 
diffuse So — > S2 photoabsorption spectrum is presented in the next section below. 

In terms of \a) states, the S2 population is 

PsAt) = §K(t)\ 2 \(a\^)\ 2 (a\Q\a), (16) 

where r a (t) = exp(-iE a t/h) sm(A a t/(2h))/(A a t/(2h)), so that |r Q (t)| 2 = [sm(A a t/(2h))/(A a t/(2h))] 2 , which tends 
to be w 1, if \ A a t/(2Ti)\ <C 1, or \t\ <C 2h/A a , which is the timescale of accuracy of coarse-grained description. So, for 
times \t\ small enough, the CW-generated S2 population is time-independent, but \a) -dependent, and can be written 
as 

P S2 = $\(a\^)\ 2 (a\Q\a), (17) 

which is fully analogous to Eq. ^ for the I7) states. 

Comparison of Eqs. (fTT]) and (fT"2|) provides the following relationship between the CW-generated Ps 2 and a(E a ) 
using |a) states: 

P ^ = <r(E a ) (a\Q\a), (18) 

which is similar to Eq. (| 1 1 1) in terms of I7) states. The term (a\Q\a) = J2 K \( K \&)\ 2 in Eq. (fT5|) is the incoherent sum, 
and is relatively slowly varying function of time. When cr(E a ) displays ORIT at specific energy E = E a due to the 
destructive interference caused by the overlapping resonances, Ps 2 is also expected to be small. This expected Ps 2 
behavior is fully confirmed computationally, using our pyrazine dynamics and control software from Ref. 



16(. This 



serves as a "theoretical high resolution spectroscopy" , justifying the importance of the ORIT in the cross section for 
the computed S2 population. As an example, the computed Ps 2 , corresponding to the total cross section a(E) in Fig. 
[2]in the vicinity of E = E™g 4 , is shown in Fig. [6] The Ps 2 profile resembles the a(E) profile in Fig. [2j the distortion 
of Ps 2 shape from <j{E) shape is mostly due to the varying (a\Q\a) factor in Eq. (fl8|) . 

Beyond theoretical interest, ORIT can be of practical importance for certain applications. For example, as in 
"EIT Spectroscopy" 22 , one can selectively excite B in an A + B molecular mixture by a CW laser while leaving A 
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FIG. 6: CW-generated 52 population Ps 2 , corresponding to the total cross section o(E) in Fig. [2] Position of -E^g 4 is shown 
as a vertical dash-dotted line and an arrow on the Energy axis. 

in the ground state. To do so, one can use the local optical transparency of A at a certain laser frequency, when A 
does not absorb, provided that B absorbs well at this laser frequency. 



V. COMPARISON OF THE COMPUTED S -»■ S 2 PHOTOABSORPTION SPECTRUM WITH THE 

EXPERIMENTAL DATA 



Finally, we estimate the accuracy of our pyrazine vibronic structure and associated approximations by computing 



the Sq — >■ S*2 photoabsorption spectrum 1(E) and comparing it with the available experimental data [21 
1(E), one can convolute the stick spectrum, composed of |(a|/i|(7)| 2 values, with Lorentzian (as in Ref. 
uniform FWrlM^ in energy domain: 

r/2 



To obtain 
23[), having 



I(E)cxu Etg J2\{a\ 



Ml: 



(19) 



(^ 9 - WQ , g ) 2 + r 2 /4' 

where i^E.g = (E — E g )/Ti, and F" = FWHM^/ft. Alternatively, 1(E) can be computed using the Fourier trans- 
form of the autocorrelation function C(t) = (*(0)|#(i)), where |*(0)) = fi\g) = |* e ), |#(t)) = U(t)\V(0)) = 
y e xp(— iEyt /fi) I7) (7I v\.q) , thus giving C(t) « ^ Q \(a\fi\g)\ 2 T a (t) in terms of coarse-grained states \a). Then, as in 



Refs. 



20 



24l | . 1(E) is obtained as 



P 00 

1(E) oc cj^gRe / dt exp(— t/T d ) C(t) exp(ioj E , g t), 
Jo 

where Tj, is the damping parameter, giving the spectral broadening. 



(20) 



11 



The spectrum 1(E), computed using Eq. ([19]) with FWHM E of 0.05 eV, and 1(E), computed using Eq. ([20 



with Td — 25.0 fs, are shown in Fig. [7] along with the experimental spectrum of Yamazaki et al. 



21|. Fig. [7] shows 



that the overall shape of the computed spectra (which are almost identical to each other) is in very good qualitative 



agreement with experiment (e.g., better than in Refs. 



23 



25j). Our computed spectra does lack the third peak near 



the maximum and somewhat overestimates absorption at low energies (before maximum), being less accurate than in 



Refs. [20J, [24| . These quantitative differences likely reflect our use of the approximate coarse-grained \a) states, the 
simplifying restriction we impose [l^] to a total of 176 S2 bright states, and the use of a phenomenological T (or T^) 
to describe the broadening of the S2 states. These approximations, however, in no way affect the general physics of 
ORIT or of constructive spectral enhancement, which are the central focus of this paper. 
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FIG. 7: Experimental pyrazine So — > S3 photoabsorption spectrum (full line) from Ref. 
spectra (long dashed line and short dashed line). Here, FWHMb = 0.05 eV, Td = 25.0 fs. 



211 ] together with two calculated 



VI. SUMMARY AND CONCLUSION 



In conclusion, we have computationally demonstrated interference-induced transparency (and cross section 
enhancement) in a polyatomic molecule, using the Sq — > S2/S1 photoabsorption cross section in pyrazine as an 
example. Specifically, transparency has been shown to arise from destructive quantum interferences associated with 
overlapping resonances [9J which are optically accessible in the So — > S2 photoexcitation. In the specific case examined 
as an example, the destructive interference between several locally important resonances was found to be responsible 
for the local transparency effect. Interference induced cross section enhancement was also shown. These effects are 
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expected to be ubiquitous in polyatomic molecules, a feature which has already been shown relevant to the coherent 
control of internal conversion 13j, [16[ . 
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Appendix: Consideration of Pyrazine Rotational Spectrum 

Pyrazine is a rigid aromatic system. For low rotational angular momentum quantum numbers, one can neglect 
rovibrational interaction and consider the molecular Hamiltonian as a sum of rotational and vibronic parts, being inde- 
pendent of one another. This gives the molecular wavefunction as a product of rotational and vibronic wavef unctions, 
and the molecular energy as a sum of rotational and vibronic energies governed by different quantum numbers. 

arazine is a nearly symmetric oblate top with rotational constants B = (A + B)/2 m 0.2 cm -1 , and C ps 0.1 
,Q- Its rotational energy can be approximated as E R (J,K) = BJ(J+1) + (C-B)K 2 ps 0.2 J(J+ 1) - 0.1K 2 
cm -1 , where J and K are rotational angular momentum and its projection on the figure axis of the top, respectively 



(J > K). The symmetric top rotational selection rules 



allow only | J", K") — > | J', K') transitions, belonging to P, 



Q, and R branches, having A J = J' — J" = — 1, 0, +1, respectively, with AK = K' — K" = 0, ±1. In case of pyrazine, 
belonging to the point symmetry group D2h, Sq{ 1 A 9 ) — » S% ( -8211) transition dipole moment belongs to B2 U symmetry 
(denoted T y or y in the B>ih character table). This transition dipole moment lies in molecular plane, perpendicular to 



the figure axis of the top 



29 



30] . thus retaining only AK = ±1 transitions in the So — > S2 rovibronic spectrum. 

If the ground state is rotationally very cold, being comprised of only one lowest rotational state | J" = 0, K" = 0) , 
then only the R branch exists. In the current case this gives only two R branch transitions \J" = 0, K" = 0) —> | J' — 
l,K' = ±1). Both transitions have relative rotational line strengths of 1/2 [28j, and the same AEr = Er(J' = 
\,K' = ±1) - E R {J" = 0, K" = 0) pa 0.3 cm" 1 , due to rotational if-doubling. Both transitions merge into the single 
line in the spectrum, uniformly shifting the vibronic spectrum by 0.3 cm -1 = 3.72 x 10~ 5 eV, which is very small in 
comparison with the energy of vibronic transition, starting here from 4.06 eV, Sq — > Si vertical electronic excitation 
energy [3l|. Considering the rotationally very cold ground state, one indeed can focus only on vibronic transitions. 
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